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Here, we study the classical information capacity of a quantum channel, assuming linear optical 
encoding, as a function of available photons and optical modes. We present a formula for general 
channel capacity and show that this capacity is achieved without requiring the use of entangling 
operations typically required for scalable universal quantum computation, e.g. KLM measurement- 
assisted transformations. As an example, we provide an explicit encoding scheme using the resources 
required of standard dense coding using two dual-rail qubits (2 photons in 4 modes). In this case, our 
protocol encodes one additional bit of information. Greater gains are expected for larger systems. 
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I. INTRODUCTION 

Superdense coding is an elegant application of basic 
quantum mechanics that can provide exciting gains in 
the capacity of an information channel. A simple pro¬ 
tocol proposed in 1992 by Bennett and Wiesner proved 
that two classical bits of information could be sent over 
a quantum channel via a single qubit [1], a compression 
made possible by exploiting entanglement as a physical 
resource. Entanglement can manifest in many forms and 
in a variety of physical hosts; engineering an entangled 
system to be used as an efficient communication chan¬ 
nel thus presents an interesting challenge. In general, 
an optimal quantum channel will require a well-devised 
protocol to properly manipulate a maximally entangled 
state [2]- 

We focus here on a linear optical [3] quantum com¬ 
munication channel constructed from a set of modes and 
photon quanta. The Hilbert space dimension for such 
a system of N total photons propagating in M optical 
modes is given by 


_ (A^ + M-l)! 
N\{M-iy. 


( 1 ) 


or the number of ways to place N photons in M modes. 
Considering, for example, a system of two dual-rail qubits 
{N = 2, M = 4), we note that dn = 10 whereas the 
logical qubit space is only four-dimensional. Resource 
efficiency is paramount in communication and such over¬ 
head immediately limits channel capacity. In this paper, 
we aim to develop optimal dense coding protocols specif¬ 
ically for a linear optical quantum channel, without im¬ 
posing the qubit as an operational basis. 

The protocol requires that Alice and Bob initially share 
an entangled state of N photons, l/’i), with Ma out of 
the total M modes under Alice’s control. Alice chooses 
a classical symbol x and encodes it by operating on her 
modes, taking |^i) to I'/a:). Her resources are forwarded 
to Bob, who then performs a quantum measurement on 
the entire state and reads out a classical bit string. Re¬ 
fer to Fig. 0. dn serves as an absolute upper bound 
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FIG. 1: (Color online) A typical quantum optical encoding 
device consisting of 4 modes. 2 modes have been distributed 
to Alice (top), and 2 modes have been distributed to Bob 
(bottom). Alice applies a linear optical operation and sends 
her resources to Bob who performs a measurement on the 
total state. 


on the number of distinguishable states Alice can send, 
though due to the restriction of utilizing only linear opti¬ 
cal operations (beam splitters and phase shifters), Alice 
may not be able to make use of the entire Hilbert space, 
even if she possesses all of the modes {Ma = M). The 
crux of the problem is that not all quantum logic oper¬ 
ations are guaranteed to be accessible via linear optical 
transformations. Entangling operations in particular re¬ 
quire photons to interact, a task they are notoriously ill- 
suited for. Solutions to this problem exist, e.g. the use 
of KLM [1] auxiliary resources or non-linear Kerr media. 
For now, we omit such components from our system. As 
we shall determine, linear optics is sufficient to construct 
a maximally efficient communication channel. 

We represent the state of our system in the Fock basis 
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A general linear operation can then be expressed as 
Ma 

d- ^ UijO^j V i = 1... Ma , (3) 

i=i 

where Uij are the matrix elements of a unitary operation 
U acting on Ma modes. Alice encodes one of |A| classical 
symbols into the system by generating the appropriate 
quantum state, , |'02), • ■ •, |'0|x|)- Each state I'l/'x) 
is constructed by applying a different transformation 
to the initial entangled state, |'0i) (here without loss of 
generality we will take Ui = /). 

The mutual entropy H{X ■. Y) shared between Alice 
and Bob is constrained by the Holevo bound, which for 
pure states can be written as 

H{X : r) <1 S{p) <2 H[X) <3 log2 |A|, (4) 

where X and Y are the encoding and decoding variables, 
1^1 

P = ^Px\i’x) (5) 

x—1 

Px is the probability of choosing symbol a;, S{p) is the 
von Neumann entropy of p, and H{X) is the Shannon 
entropy of Alice’s encoding variable. Inequality 1 reflects 
the fact that depending on Bob’s measurement capabil¬ 
ities, he may not be able to access all the information 
encoded by Alice; it becomes an equality only if Bob can 
perform arbitrary measurements on words of arbitrary 
length [SHE]- We return to this point in Sec. VI, and in 
the following we focus on S'(p), or the encoding capacity 
of the channel. Inequality 2 becomes an equality when 
all contributions \'4>x) to p are orthogonal, and finally in¬ 
equality 3 becomes an equality when all symbols in X 
are chosen with equal probability. 


II. AN UPPER BOUND ON ENCODING 
CAPACITY 

From Eq. Q , it immediately follows that encoding ca¬ 
pacity is limited by the Hilbert space dimension, 

S{p) < log2(dff). (6) 

This inequality provides a valuable reference point, but 
for our purposes we can derive an even stricter bound on 
S{p). First, we abstractly define ds as the dimension of 
the span of the states Alice can generate via linear opti¬ 
cal transformations as presented in Eq. (|^. The Holevo 
bound again restricts the encoding capacity, 

S{p) < log2(ds) ds <dH ■ (7) 

Then, we note the Hilbert space of our channel can be 
decomposed into a direct sum 

N 

•H = 0 (8) 

Na=0 


where each term represents the subspace corresponding 
to a fixed distribution of photons between Alice and Bob. 
Na is the number of photons in Alice’s modes, and Nb = 
N — Na- Any initial state \tpi) can be decomposed into 
a sum of components 

IV'jva) e . (9) 


For each such component we can construct a Schmidt 
decomposition over a set of min(d^’^'^, ) basis vec¬ 
tors, where d^’^^ and are the dimensions of 

'^a,Na and respectively. An?/operation on Alice’s 

modes will act only on the basis vectors of her subspace. 
It follows that for any where < d^’'^®, Al¬ 

ice has no control over the dimensions of Bob’s subspace 
which are not spanned by the basis vectors of the Schmidt 
decomposition. Thus, the upper bound for the dimension 
of the span of the set of states Alice is able to generate 
within each term of Eq. (|^ is given by 

d^^ < d^^'^ min(d^’^'", d®’^^). (10) 


We combine Eqs. 


the total span ds for given p 
Ma,Mb > 1 : 


and (10) to find the upper bound on 
loton number N and modes 


ds < 


N 


E 

Na=0 


( 11 ) 


where 


f{NA) = g{NA,MA)mm{g{NA,MA),g{NB,MB)), ( 12 ) 


g{n,m) 


(n + m — 1)! 
nl{m — I)! 


(13) 


is the number of ways to distribute n photons over m 
modes, and Mb = M — Ma- 

This is a key result; if we combine Eqs. Q and 0 
we hnd an upper bound on the information capacity for 
a given physical device. As we will see in Sec. HI, we 
find that the bound obtained here is tight. 


III. NUMERICAL TESTING OE ALICE’S 
ENCODING CAPABILITIES 

Eqs. (|7j) and 0 provide a fully analytic upper bound 
for the information capacity of an optical channel, but the 
question remains whether and how this bound can actu¬ 
ally be reached using linear optical encoding. To obtain 
the actual encoding capacity, we numerically maximize 
S{p) over the initial state \ijji) and Alice’s transforma¬ 
tions Ux- We call the result of this optimization S-^ax- 
Similarly, ds is computed by evaluating the rank of the 
density matrix for an ensemble of states generated with 
a sufficiently large set of random Ux- 

For small N and M, the numerical results are fairly 
straightforward and elegant. Of particular interest is the 
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case where TV = 2, M = 4, Ma = 2, which uses the 
same physical resources as the standard, two dual-rail 
qubit dense coding protocol. Here we find that ds = 8, 
in agreement with the right hand side of Eq. ( |TT| ) . Fur¬ 
thermore, we find that Smax —t log 2 (ds) = logJ(8) bits 
can be achieved by generating |X| = 8 completely dis¬ 
tinguishable states, see Fig. (Left). In contrast with 
the standard dense coding protocol, which allows Alice 
to communicate 2 classical bits via a single photon in 
two modes, we instead find that Alice is able to send 
3 classical bits via 1.25 photons (2.4 bits per photon or 
1.5 bits per optical mode). This is an important result; 
Alice is able to encode additional information into the 
channel by not restricting herself to the qubit basis. The 
advantage gained by expanding the optimization space 
beyond the dual rail basis has been consistently observed 
in other work on gate optimization and photonic quan¬ 
tum communication |10j . 
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Instead, she must generate a set of at least 20 non- 
orthogonal states to obtain the maximum capacity of the 
channel. It becomes apparent that Alice cannot, in gen¬ 
eral, construct a set of ds orthonormal basis vectors over 
the entire subspace accessible to her via allowed opti¬ 
cal transformations. She must instead saturate the sys¬ 
tem with additional linearly independent states until she 
reaches the potential of the hardware. 

In Fig. [§ we present a larger system, N = 3, M = 6, 
Ma = 3, where ds = 38. Here the full encoding capacity 
log 2 ds is approached slowly as the number of symbols 
increases, although a modest number of symbols is suf¬ 
ficient to attain a very large fraction of the maximum 
capacity. 
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FIG. 2: (Color online) Maximization of the von Neumann en¬ 
tropy S{p) as a function of the number of symbols |A|. (Left) 
N — 2, M = A, Ma = 2. A clear convergence to log 2 ( 8 ) = 3 
is observed. (Right) A = 3, M = 5, Ma = 2. Here Smax 
reaches log 2 |A| for |A| < 13, indicating that up to 13 or¬ 
thogonal states may be generated. Globally, Smax reaches 
log2(18) at |A| = 20. Analytic solutions can be constructed 
for |A| < 13. 


For larger, more resource-intensive systems, such 
straightforward convergence of Smax to log 2 |^| is not 
observed. We instead find an intermediate optimization 
regime where the encoding capacity is achieved by in¬ 
creasing the number of coded states |A| beyond the num¬ 
ber of orthogonal states Alice can send. As an example, 
we present the results for TV = 3, M = 5, Ma = 2 in 
Fig. 1^ (Right). Here, we find dg = 18 [again in agree¬ 
ment with the right hand side of Eq. (11)] and indeed 
Smax log 2 (ds) = log2(18), but Only when |A| > 20. 
The maximum channel capacity is equivalent to one asso¬ 
ciated with 18 orthogonal quantum states, however Al¬ 
ice cannot generate 18 such states using linear optics. 


FIG. 3: (Golor online) Maximization of the von Neumann 
entropy for N — 3, M = 6, Ma ~ 3. \X\ must be increased 
well beyond ds = 38 for the von Neumann entropy Smax 
to maximize to log2(38). No gain in information capacity 
is observed if we allow Alice to send symbols from X with 
non-uniform probability. 


Empirically, we observe in all cases that ds and Smax 
do attain the upper bounds defined in Eq. (11) and 
Eq. 0, respectively, assuming only linear optical encod¬ 
ing. We therefore define classical information capacity C 
in bits per channel use as a function of TV, TVf, Ma'- 


N 


C{N,M,Ma) = log 2 (cis) = log 2 5] /{Na) , (14) 


\Na=0 


where /{Na) is defined in Eq. (12). Because the bound 
in Eq. 0 holds for any photon-number preserving op¬ 
eration on Alice’s modes, the use of non-deterministic 
or non-linear entangling optical components will not im¬ 
prove channel capacity. 
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IV. ASYMPTOTIC BEHAVIOR OF THE 
ENCODING CAPACITY 


To gain some insight into the parameter dependence 
implied by Eq. (14), we consider the thermodynamic 
limit M = aN, for Ma,Na_,Mb,Nb,M,N ^ 1. Ap¬ 
plying Stirling’s approximation, the full Hilbert space di¬ 
mension is given by 


logsdff ^ A^[log2(l + a) + alog2(l + a ^)] ■ (15) 

Similarly, 

77 - 777 

^og2 g{ri,m) ~ TOlog2(l H-) -Hnlog2(l H-), (16) 
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FIG. 4: (Color online) Encoding capacity log 2 ds for M — 2N 
and different ratios Ma/Mb- As an upper and lower bound, 
we include, respectively, the Hilbert space dimension and the 
capacity of traditional dense coding using two dual-rail qubits 
(which assumes Ma = Mb)- 


and thus 7V^, the maximum of g{NA, MA)g{NB, Mb), is 
given by 

V. DEVELOPING EXPLICIT ENCODING 
PROTOCOLS 


N^-NMa/M. (17) 


Now /(Na) is piecewise defined with the crossover point 
Na determined by Ma- One can verify that if 

^ > 1 , = VI if ^ = 1 , and < VI if ^ < 

1. Thus, the value of f{NA)max is dependent on the 
ratio of modes Ma/Mb- We find the following regimes 
of behavior as a function of the mode ratio Ma/Mb- 


> I ^ > Na ^ ds ^ dn (18) 

^ = 1^V-Vl = f (19) 

^ < 1 ^ V^ < VI ds « dff (20) 

Mb 

MB>g{N-l,MA) 

ds ^ g^{N — 1, Ma) dn (21) 


In particular, we see from Eq. (18) that Alice main¬ 


tains near absolute control of the channel in the range 
Ma/Mb > 1 , and the maximum encoding capacity is 
asymptotically indistinguishable from the capacity im¬ 
plied by the total Hilbert space dimension dn- Eq. (19) 


indicates that when the modes are evenly split be¬ 
tween Alice and Bob, the maximum encoding capacity 
is asymptotically lower by precisely one bit than the ca¬ 
pacity implied by ds- From Eq. , we find that for 
a fixed number of modes Ma controlled by Alice, any 
modes Bob possesses in excess of g{N — 1,Ma) do not 
contribute to channel capacity. 


Thus far we have presented a general method for cal¬ 
culating the information capacity of linear optical hard¬ 
ware. An important task remains; we must be able to 
design a dense coding protocol for actual implementa¬ 
tion. We return to the prototypical case: V = 2, M = 4, 
Ma = 2. Using the results of the numerical optimization 
of S{p), we find the general set of initial states and trans¬ 
formations that satisfy {ijji\tpj) = OV i j- The result is 
a non-trivial set of equations, which simplify elegantly if 
we slightly restrict the space of Alice’s transformations, 
t/ 3 , C/4, ..., Us, see Eq. ((^. 


VI. DISCUSSION 

We have already noted that the Holevo bound is not 
automatically reachable, and the information capacity of 
an optical quantum channel is limited by factors other 
than efficient encoding. The susceptibility of entangled 
states to decoherence and noise must be considered. Op¬ 
tical quantum computers are an attractive idea in this re¬ 
gard because photons do not strongly interact with most 
matter; photonic states can thus be carried over large 
distances nm. 

So far we have not addressed the practicality of 
generating the pre-distributed maximally entangled ini¬ 
tial states. The question remains as to whether 
can be constructed using deterministic optical compo¬ 
nents. This will likely depend on the physical resources 
{N, M, Ma) of a particular device. While deterministic 
state generation is certainly preferred, it is not essential 
for efficient communication if one considers cost as that 
associated with channel usage. 
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At the receiving end of a quantum channel, signal de¬ 
coding is an especially challenging task. For the case in 
which Alice sends orthogonal states, the Holevo bound is 
achieved by applying an appropriate von Neumann mea¬ 
surement. In principle, a sufficiently large subset of von 
Neumann measurements can be carried out by having 
Bob rotate incoming states into the photo-counting ba¬ 
sis states (see Lougovski and Uskov m)- We hnd that 
it is not possible to implement such a measurement us¬ 
ing deterministic linear optics exclusively. If one were to 
impose a linear optical measurement. Bob would be un¬ 
able to fully distinguish incoming states m- To make 
things even more complicated, our results suggest that 
in order to optimize S{p), Alice will often send a set of 
non-orthogonal states. Thus, even if we were able to ap¬ 


ply an optimal POVM, H{X : Y) would fall well below 
S (p) for a single communication O [12] . Achieving the 
Holevo bound for an a priori set of non-distinguishable 
quantum states is a common challenge for quantum infor¬ 
mation theory, and has been the subject of some study. 
Proposed solutions [BHHl [T3| rely on sending code words 
comprised of multiple states, For example, the 

PGM [5] scheme guarantees that a word composed of I 
states carries H{X : Y) = lS{p) information for a coding 
vocabulary in the limit of large 1. The caveat is that Bob 
must store previously sent letters in a quantum mem¬ 
ory storage device, so that he may act on an entire word 
with a joint measurement. This could be implemented 
using IM storage modes and a rail-switching junction for 
incoming letters. 


Initial entangled state: 

|2000) -k |II00) -f |10I0) -k 046 *'^" |I001) -k 056 *^'= |0200) 
|0I10) -k 076*^'" |010I) -k cge*''® |0020) -k |001I) -k cioe*‘^i« |0002) . 


Alice s operations: 


Restricted by: 




= 



-\/fe - 1 

(x = 3,4, ..., 8 ) 

= <73 + - 3) (a: = 4, ..., 8 ) 


c? 




, 3 

: = g C3 


2,2 ^ 2 , 2,2 1 

^ — C7 C4 — Ce <^3 + <^4 — ^ <^8 + <^9 + ^10 — g 

ds + dr — d 4 — de = mir m € 'Lodd 

C1C2 cos(di - d2- <73) - CiCr, sin(di - ds - 293) - C2C5 cos(d2 - ds - 93) = 

C1C2 sin(di - dg - 93) - c:i<^5 cos(di - da - 2^3) - C2C5 sin(d2 - da - q^) = 


( 22 ) 


VII. CONCLUSION 


In summary, we have shown that hardware-specific 
coding algorithms must be developed in order to make 
full use of a linear optical quantum channel. In Eq. (141, 


we obtained a simple expression for channel capacity as 
a function of photons and optical modes. Furthermore, 
we found this capacity is achievable via some encoding 
procedure using linear optics exclusively. We observe 
that for small systems, Alice can send completely distin¬ 
guishable states. Generally, however, Alice must choose 
from a pool of non-orthogonal states. In either case, 
it is possible to extract an encoding protocol from the 
numerical maximization of the von Neumann entropy. 
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